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THE  FAE  FIELD  OF  A  ROCKET  EXHAUST  JET 


AT  LOW  AND  MODERATE  ALTITUDES 

F.  P.  Boynton 

ABSTRACT 

A  method  is  given  for  predicting  the  flow  field  of  a  rocket  exhaust 
plume  at  distances  far  removed  from  the  nozzle  exit  at  altitudes  where 
afterburning  of  the  exhaust  is  appreciable.  The  calculation  combines  a 
fluid  mechanical  analysis  of  turbulent  mixing  due  to  Libby  with  the 
adiabatic  flame  temperature  calculation  of  Boynton  and  Neu.  The  eddy 
viscosity  is  determined  from  a  consideration  of  conditions  under  which 
compressible  turbulent  flows  appear  to  exhibit  self -preserving  behavior. 
Instructions  are  given  for  preparing  input  to  two  computer  programs  which 
are  based  on  the  analysis  in  this  report.  In  two  appendices,  it  is  shown 
how  the  eddy  viscosity  constants  may  be  derived  from  incompressible  jet 
flow  data  and  a  comparison  of  the  results  of  the  present  calculation  with 
experimental  wind-tunnel  rocket  exhaust  behavior  is  given. 


V 


INTRODUCTION 


The  flow  far  downstream  from  the  nozzle  exit  of  a  rocket  engine  is  governed 
primarily  by  viscous  (turbulent)  and  chemical  interactions  with  the  surrounding 
atmosphere.  The  inviscid  interactions  caused  by  an  initial  imbalance  between 
nozzle  exit  and  ambient  pressures  have  largely  been  dissipated  in  this  region,  and 
it  is  therefore  appropriate  to  consider  that  these  mixing  and  burning  phenomena 
occur  at  relatively  constant  pressure.  A  description  of  this  flow  is  then  given  by 
solving  an  Inltlal-value  problem  for  the  boundary  layer  equations .  The  initial 
conditions  are  those  of  the  rocket  exhaust  expanded  to  ambient  pressure .  For 
nozzle  exit  pressures  close  to  ambient,  the  properties  of  this  equivalent  jet  are 
given  by  an  isentropic  expansion  of  the  exhaust .  For  nozzle  exit  pressures  much 
greater  than  ambient,  the  non-isentropic  analysis  of  Thomson^  may  be  applied.  For 
nozzle  exit  pressures  much  less  than  ambient  (which  are  perhaps  of  lesser  interest) 
the  initial  properties  might  be  taken  as  those  behind  a  normal  shock  at  the  jet 
exit. 

In  the  following  paragraphs  we  shall  discuss  briefly  some  elementary  consider¬ 
ations  of  the  structure  of  the  far  field  which  may  be  derived  from  a  similarity 
analysis  of  jet  flows  in  a  coordinate  system  in  which  the  boundary  layer  equation 

reduces  to  their  incompressible  form.  We  shall  then  discuss  a  linearized  solution 
2 

obtained  by  Libby,  and  show  how  this  analysis,  combined  with  the  considerations 
obtained  from  the  similarity  analysis,  may  be  used  to  treat  problems  of  rocket 
exhaust  jets . 
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1 .  SIMILARITY  CONSIDERATIONS 

It  is  well  known  that  the  flow  of  a  turbulent  Jet  of  constant -density  fluid 
far  from  the  origin  is  self-preserving :  the  variation  of  any  mean  quantity  over 
any  plane,  x  =  constant,  is  expressible  non-dimensionally  through  suitable  scales 
of  length  and  velocity,  and  u^,  as  a  universal  function  of  Analyses  of 

the  equations  of  motion  for  either  the  case  where  the  jet  velocity  is  much  greater  . 
than  that  of  the  ambient  medium  or  the  case  where  the  Jet  and  ambient  velocities 
differ  by  a  small  amount  have  shown  that  this  behavior  is  possible,  and  a  large 
number  of  experiments  bear  out  that  it  really  occurs.^  However,  the  existence  of 
full  self-preservation  in  a  Jet  flow  with  large  density  variations  leads  to  an 
absurd  result,  that  the  density  varies  as  a  power  of  the  axial  distance.  While 
this  result  is  possible  over  a  small  region  of  the  flow,  its  continuance  violates 
the  condition  that  the  density  must  approach  its  free-stream  value  far  downstream. 

A  transformation  of  the  boundary  layer  equations  exists  which  allows  the 
coordinate  system  to  be  stretched  or  squeezed  in  such  a  manner  as  to  eliminate  the 
density  as  a  variable.  The  equations  then  assume  their  incompressible  form,  and  it 
is  possible  for  self -preserving  profiles  to  exist  in  the  new  coordinate  system.  The 
transformation  is  that  of  Howarth^  and  Dorodnitsyn,^  wherein  a  new  radial  coordinate 
is  introduced  as 


R  = 


o 


(1) 


where  is  a  suitably  chosen  reference  value  of  the  density.  The  boundary  layer 
equations  in  terms  of  mean  flow  properties  are  these  (for  axially  symmetric  flow): 


2 


9u 


Continuity:  ^  (rpu)  +  ^  (rpv)  =  0 


Momentum:  pu|^  +  [j.  T  , 

ox  ^  dr  5x  r  dr 

2 


dr  5x  ^  r  dr  '•*  xr" 

2 

1  =  i  ^ 

I’  '  '*r  xr' 


Energy;  pu  g-Che?^]  +  py  ,  i  ^  {r(9  V^)) 

Species  continuity:  pu  ^  t  pv  ^  .  i  ^ 

The  transfonnation  (x^r)-.(x,R)  is  performed  by  noting  that 

dx  dx  dx  dR^  dr  “  dr  ^ 


and  by  requiring  the  stream  function  Y  to  be  the 


same  in  each  flow , 


* 

Y  =  Y 


By  using  the  relations  (3)^  and  also  noting  that 


(2) 


(3) 


*  1  dY  * 

U  =  —  ^rr.  V  = 


R  dR' 


1  ^ 
R  dx" 


(5) 


whence 


*  r  pv  dR 

R  Po  dx 


(6) 


the  following  equations  are  obtained  in  place  of  the  boundary  layer  equations : 

Continuity:  JJ  +  i  ^  (Ey*)  =,  0 

ttaentm.:  u  |S  ,  /  |i  ,  ^  ^ 

Energy:  u  ^Chtif]  t  v*  „yT*)jl  (7) 

dc.  up  dc. 

Species  continuity:  u  -rri  ^  ■y*  _ i  _  1 

dx  ^  ^  dR  -  p  R  dR^^^ri-® 
o 

The  reader  will  recognize  that  these  are  simply  the  boundary  layer  equations  for  an 
equivalent  incompressible  flow.  The  shear  forces  and  flows  of  heat  and  matter  are 


3 


required  to  be  the  same  for  unit  mass  in  both  flows,  whence 


TO  * 

R  T  p  dV 
xr  o 


TO  * 

R  q  p  dV 
^r  o 


TO  * 

R  J  .  P  dV 


T 

r  T  p  dV 
xr 


T 

r  q  p  dV 

r 


T 

r  J  .  P  dV 


where  dV  and  dV  are  volume  elements  in  the  transformed  and  original  coordinate 
systems . 

We  now  wish  to  express  the  set  of  equations  (7)  in  a  self -preserving  form.  We 
can  do  this  in  two  limiting  cases,  just  as  for  incompressible  flow.  If  we  let  u^ 
be  the  velocity  of  the  external  stream,  and  let  u^^  be  the  difference  between  the 
local  velocity  in  the  jet  and  the  free  stream  velocity,  we  my  write  the  momentum 
equation  as 


(u  +  u^  ) 
'  e  1' 


^^1  *  ^^1 
dx  ^  ^  dR 


^  A  fT? 
p  R  dR  ^  '^'xr'^ 


(9) 


We  can  either  assume  that  u^  »  u  at  every  location  of  interest  to  us,  or  that 

1  e 

u  »  u, .  These  two  cases  will  be  investigated  separately, 
e  1 


Case  A  (u^  »  u  ) : 

X  6 

Here  we  choose  the  self -preserving  functions  u^  =  u^  ^(^1),  v  =  C(Tl)  and 

A  =  p  u  ^  g(T]),  where  Tl  =  B./&  .  We  substitute  these  relations  into  equation  (9) 
xr  o  m  '  m 

and  perform  a  second  coordinate  transformation,  (x,r)~*(x,1l) .  The  result  is 

2 

il 


du 


u 


m  dx 


“A 


u  ^  d^  u 

m  m  _  m 


i  dx 
m 


1. 


m 


^  (%)'  - 


(10) 


The  requirement  for  self-preservation  is  that  the  coefficients  of  the  universal 


functions  of  T]  be  either  zero  or  proportional  to  each  other.  This  requirement  is, 

1  “S 

satisfied  if  £  =  -  (x-x  )  and  u  ~  (x-x  )  ,  where  a  is  at  most  a  function  of  the 

m  O'  o  m  o 


h 


initial  conditions  and  a  is  not  yet  determined.  The  momentum  Integral  conservation 
equation  is 

oo  CO  CO 

2n  rdr  =  2p  tt  1  RdE  =  ttu^  p  [2  f^TldTj!  ^  nr^p.u.  ,  (ll) 

1  oJl  mm  oLJ  J  j  ^  ^  ' 

o  o  o 

which  is  obtained  by  letting  the  momentum  flux  across  an  axial  station  be  that 
through  the  nozzle  exit.  The  bracketed  tenn  is  a  constant,  and  we  can  set 


This  is  the  same  relation  (except  for  the  presence  of  the  term  containing  the  nozzle 
and  reference  densities)  as  is  obtained  by  Hinze  for  an  incompressible  flow.  Thus 
a  =  1,  and  u  ~  (x-x  ) 

The  form  of  equation  (l2)  suggests  that  the  universal  function  f  (Tl)  may  be  the 
same  for  compressible  and  incompressible  jets .  Data  from  several  sources  have  been 
analyzed^  to  determine  whether  this  congruency  exists .  It  was  found  that  for  jets 
in  which  the  >fe,ch  number  was  low,  but  where  density  variations  were  large,  the 
profiles  fell  on  the  same  curve  when  the  reference  density  was  chosen  to  be  p^ . 
Efata  from  jets  in  which  the  initial  Jfech  number  was  high  gave  curves  which  were 
narrower  than  the  low-speed  curve,  but  which  could  be  brought  into  congruence  with 
the  low-speed  curve  by  assuming  that  ct  was  a  function  of  the  initial  tfech  number. 

A  choice  which  fit  the  data  examined  (which  included  only  initially  isoenergetlc 
jets)  quite  well  was  cr  =  /l  +  i(Y-l)M^'.  This  analysis  indicated  that  the  axial 
velocity  decay  of  all  jets  in  a  stagnant  ambient  medium  would  be  the  same  when 
plotted  in  an  axial  coordinate  proportional  to 

”0  "b  ’ 

and  all  jets  examined  followed  this  behavior.  (We  might  note  that  this  behavior 


5 


Y 

was  hypothesized  by  Thring  and  Newby'  in  a  study  of  turbulent  flames.) 

To  this  point  we  have  said  nothing  about  the  shape  of  the  self -preserving 
functions  f(1l),  §(11)  and  g(Tl)  except  that  they  can  be  made  the  same  for  all  Jets. 

In  order  to  proceed  further,  we  require  an  assumption  relating  the  shear  stress  to 
other  parameters.  Such  a  relation  is  obtained  by  introducing  the  concept  of  an 
eddy  viscosity,  e^,  which  is  defined  such  that 

.  du  T  .  . 

5?  '  V  d3) 

In  many  incompressible  shear  flows  reasonable  agreement  of  calculated  and  experi¬ 
mental  velocity  profiles  are  obtained  by  assuming  is  constant  across  any  plane 
of  constant  x  in  the  flow.  The  congruence  of  incompressible  and  (transformed) 
compressible  velocity  profiles  suggests  that  a  similar  assumption  could  be  applied 

to  the  compressible  flow  if  the  value  of  e*  for  the  transformed  coordinates  is  used. 

* 

The  -values  of  and  may  be  related  to  each  other  by  remembering  that  the  momentum 
flow  per  unit  mass  (l.e.,  the  force  exerted  by  shear  stresses  on  a  unit  mass)  is  the 
same  in  the  two  systems .  Thus 

m 

RT^r  Pq  prdrdxdG  (l4) 


or 


Wow,  since 


and  so 


e*  p  ^  ^  dE 
V  o  3R 


2 


P 


5  dR  5 

3r  Sr  Sr  ^ 

5u  SR  Su 

Sr  “  Sr  SR  ^ 


e 


V 


* 

e 

V 


(15) 


(16) 
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This  representation  of  the  eddy  viscosity  was  first  suggested  in  print  by  Ting  and 

T  -KK  ® 

Labby . 

* 

For  actually  calculating  shear  flows^  we  require  an  expression  for  e^.  Since 


e  is  defined  as 
V 


* 

s  = 
V 


ip* 

xr 


D  ^ 

o  SR 


(17) 


T*  * 

we  may  substitute  the  self -preserving  function  for  T  and  u,  and  express  e  as 

XT  ’  V 


V  m  a 


(18) 


* 


In  order  for  to  be  constant  for  a  given  value  of  x,  the  ratio  g('Tl)/f  ^  (Tj)  must  be 

a  constant,  K.  If  we  identify  u  with  u  ,  the  value  of  u  on  the  centerline,  and 

recognize  that  the  transformed  jet  half -radius,  Ri,  is  simply  T|i  (x-x  ),  we  may 

■&  s  0 

* 

express  e  as 


e  =  (k/ctT1i)  Ri  u  . 

V  e  3  ^ 


(19) 


This  is  the  form  generally  used  for  in  incompressible  jet  flows  in  a  stagnant 


medium,  with  the  exception  of  the  "Ifech  number  correction",  a. 

If  it  is  assumed  that  is  constant  across  a  plane  of  constant  x,  then 
Equation  (lO)  may  be  solved  for  f .  The  result  is 


f(''l)  =  ^ 


(20) 


* 


V 


where  o'  ~  — .  This  solution  has  been  found  to  fit  experimental  radial  velocity 
c 

profiles  in  incompressible  jets  very  well  except  in  the  outermost  regions  of  the  jet, 
where  the  calculated  velocity  drops  off  less  rapidly  than  the  measured  velocity. 


(!te,se  B  (u  »  U-,  ) : 

e  1' 

^y. 

Again  we  introduce  the  self -preserving  functions  u^  =  u^f  (T|),  V  =  u  C(T1), 


*Wote  that  R^  is  proportional  to  a  displacement  thickness . 


2 

P^gCIl).  The  momentum  equation  becomes 


pdu  u  Ai 

u  l^f  “  “  T]f' 

<=  I  dx  , 


m 

j^m 


e  Ldx 

The  integral  equation  of  momentum  conservation  is 

OO  CO  03 

,2 


(21) 


2j  PUi(ui4-u^)rdr  =  2^  p^u^  (u^+u^  )RdR  =  2p^^; 


u  u 
m  m  e, 


f  TldT]  =  u  u .  p  .  r  . 

e  j  " 


(22) 


and  the  combination  of  (2l)  and  (22)  establishes  that  for  self -preserving  flow 


1/3 

~  ^  ^  ^  •  In  this  case  we  may  express  u,  as 


-2/3 


1 - n 

U,  =r  -i  (u.-u  )  - si 

\  n  (2>'  vr_/ 


2/3 


1  -  u  “e^  iL 

^  <_  n 


(23) 


where  we  have  again  chosen  p^  as  the  reference  density. 

If  we  introduce  an  eddy  viscosity  which  is  constant  across 
then  the  equation  for  f  may  be  solved  to  show 


a  plane  x  =  constant ^ 


(24) 


Again^  it  is  possible  to  express  as 


e  =  (const.)  X  (u  Ri/o) 

V  Co 


(25) 


The  Jfech  number  correction  factor,  o,  seems  to  reflect  a  dependence  on  the 
partition  of  the  total  enthalpy  of  the  jet  into  thermal  and  kinetic  energy.  In 
these  terms,  it  might  be  more  reasonable  to  speak  of  it  as  a  Crocco  number  effect. 
The  origin  of  this  effect  is  not  clear,  but  it  may  have  to  do  with  the  production 
of  sound  by  the  turbulent  jet.  A  recent  review  by  Ribner^  quotes  results  derived 
by  Haillips,^*^  who  found  an  expression  for  the  acoustic  efficiency  of  a  mixing  layer 
at  very  high  Whch  numbers  which  may  be  expressed  as 
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71 


acoustic  energy  flux  _  w-3/2  iF 

pu3  ?  7  ‘ 

Here  2U  is  the  velocity  difference  between  the  streams ^  M  is  the  tfech  number,  v^ 
is  the  mean  square  component  of  the  fluctuating  velocity  normal  to  the  flow 
direction,  is  a  length  scale  of  these  fluctuations,  and  L  is  proportional  to 
the  thickness  of  the  mixing  zone.  Ribner  quotes  experimental  results  on  jet  noise 
from  rockets  and  aircraft  engines  which  seem  to  indicate  that  the  acoustic  efficiency 
at  high  exhaust  velocities  is  either  independent  of  jet  velocity  or  decreases  slowly 
with  increasing  jet  velocity. 

£1  order  to  maintain  a  constant  acoustic  efficiency  (in  the  context  of  Phillips' 
asymptotic  expression)  as  the  Jfech  number  is  increased,  the  ratio  of  the  relative 
turbulent  intensity  vVu  to  the  relative  scale  i/L  must  vary  as  The  eddy 

diffusivity  of  a  turbulent  medium,  on  the  other  hand,  may  be  written  as  e  ~  v^f,  or 
as  the  product  of  a  (rms)  fluctuating  velocity  and  a  scale.  The  asymptotic  expression 
for  a  suggested  by  the  jet  data  is  e  ~  M“^,  and  the  variation  of  ct  in  the  regime 

3A 

studied  is  more  like  M  .  This  is  compatible  with  Phillips  *  theory  if  v^/u  remains 
more  or  less  constant  while  f/l  decreases  with  increasing  Mich  number.  This  implies 
that  the  eddies  are  reduced  in  size,  perhaps  by  the  compressive  action  of  shock  waves, 
as  the  Mach  number  is  increased.  The  whole  field  of  turbulence  at  high  J&ch  numbers 
remains  quite  speculative,  and  considerably  more  theoretical  and  experimental  study 
is  needed  before  any  firm  conclusions  can  be  drawn.  In  the  meantime,  we  must  make 
do  with  the  skimpy  amount  of  data  available .  (The  existing  measurements  of  turbulence 
quantities  in  jets,  such  as  those  of  Laurence,  have  only  been  made  up  to  M  =  0.8.) 

2.  LIBBY'S  THEORY 

Self -preserving  flow  is  only  found  in  the  far  field  of  Jets  when  the  ambient 

9 


is  essentially  stagnant,  or  v,hen  the  difference  between  Jet  and  aabient 

velocities  is  saall.  m  order  to  develop  a  method  which  conld  be  extended  to  other 

parts  of  the  Jet  and  which  can  be  used  to  treat  flows  where  is  neither  much 

greater  than  u^  nor  nearly  equal  to  libby  obtained  a  linearised  solution  to  the 

equations  of  motion  of  a  Jet.  The  momentum  equation,  with  the  shear  stress  expressed 
by  Equation  (13),  is 


Sx 


9u  1  S 


Su> 


3r  r  3r  Sr' 


(26) 


The  procedures  for  obtaining  the  solution  are  slightly  different  when  u^  is  finite 
and  when  u^  =  0,  but  the  solution  for  finite  u^  reduces  to  that  for  u  ,0  in  the 

T  «  •  ^ 

imit  When  transformed  to  the  physical  plane.  We  shall  consider  first  Case  B,  u 
finite,  and  then  (^se  A,  =  0.  These  correspond  to  the  asymptotic  cases 
’^1  \  ■the  self -preserving  solutions. 


Case  B  (u  ^  O) : 


where 


We  apply  the  von  Mises  transformation  to  Equation  (26),  l.e.,  (x,r)-^(x,^). 


pur  =  p  u  ij; 


pvr  = 


-  p  U  (1; 

e  e^ 


9i|; 


and  obtain  the  result 


(27) 


where 


—  _  1  5 

3x  “  U^\|;  dij;  L 


ep^r^u/p/u^,^),  |Hj 


r  = 


t  P  u  I 

2  J  iv 


(28) 


(29) 


At  this  point  we  introduc 


e  the  eddy  viscosity  in  the  Howarth-Dorodnitsyn  coordinates, 
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which  may  be  shown  to  be 


*  2  2  ,  2„2 

e  z^prc/p  R  , 

V  V  O  ^ 


*  2  2 
e  =  p  r  e 
V  v 


T 

/P  P  u 
''  o  e  e  J 


2ilf^d\|f^ 


(30) 


Thus  Equation  (28)  becomes 


5u  I  5 
dx  “  u  i|f  Sil? 


1,2  J 


2ilt'dilf' 

u 


a^J 


(31) 


* 


If  we  now  assume  that  is  independent  of  r  (and  thus  of  i|r)^  and  make  the 


assumption  that  the  term  in  curly  brackets  is  identically  one^  we  find  that 


^  P  S  (X  )  "N  ^ 

ou  _  _o  V^  ‘  O  p  OU-, 

Hx  "  p  u  ■ij  dt];  df-* 


(32) 


The  assumption  regarding  the  integral  is  essentially  a  linearization,  for  the 
bracketed  term  is  really  one  only  when  u-u^  «  u^ .  We  shall  examine  the  effects  of 
the  linearization  later. 


We  now  transform  the  axial  coordinate  according  to 

X  * 

0  p 

^  V  '^o  , 

^  -J  ’ 

^  e  1  e 


(33) 


where 


P  .u .  1 
il/  .  =  r  .  (-J— i)2 


We  then  obtain  the  heat  conduction  equation. 


m 


with  the  boundary  conditions 
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u  . 

e 


u(0,ij;)  =  u 
u(0,^)  ^ 
lim  u(§,i|;)  = 


J 


The  solution  to  this  problem  is  given  by  Carslaw  and  Jaeger^  as 


u-u 


u.-u 
J  e 


-  =  (l-e 


\  *  /  1 
)  P  ( 


t/t. 


d  J 


;)  . 


(35) 


where  the  function  P*  (|  ,  |)  has  been  tabulated  by  Jfesters.^^ 

In  order  to  transform  to  the  physical  plane,  one  must  specify  the  density  as  a 

'X* 

function  of  ([f  and  e  and  perhaps  as  a  function  of  x.  It  was  found  from  the 
similarity  analysis  that  was  a  reasonable  choice,  since  it  allowed  Jets  of 

varying  density  to  be  expressed  by  the  same  radial  profiles.  The  expression  for  e* 

'  V 

may  be  obtained  by  considering  that  in  the  upstream  portion  of  the  Jet,  the  behavior 
IS  similar  to  that  of  a  two-dimensional  mixing  layer,  and  in  the  downstream  portion 
the  behavior  is  that  of  a  fully  developed  axially  symmetric  Jet  flow.  Therefore,  in 
the  upstream  region 


e 


* 

v 


Cl  x(u.-uj/a  , 


(36) 


and  in  the  downstream  region 


(37) 


The  transition  point  may  be  chosen  to  give  good  agreement  with  experiment. 
Therefore,  we  may  divide  the  transformation 


X 


(i  § 


/ 
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into  two  parts .  Above  a  certain  value  5  =  5  ,  we  have 


C_  U.-U  f 

g  =  f  -1  dx' 

^  J  CT  U  llf  . 

o  ®  J 


or 


t  .  u .-u  „ 

^  _ e  /2^\2 

2  a  u  ^ilf  '  J 
e 


r  o  21 


X  = 


IS 


Lu,- 


.  U  p  r.. 
J  e  "^e  1  j 


(36) 


Beyond  we  have 


X  + 

o 


CTu  il;  .d? 

e^J 


(39) 


This  relation  may  either  be  evaluated  directly,  or  we  may  develop  a  somewhat 
more  general  relation  which  will  allow  changing  Cg  and  <t  for  a  given  computation 
more  easily.  By  introducing  the  variable  S  =  2%/^/.,  we  can  write  (39)  as 


ilr  .  p  au  \lr  .dS^ 

J  I  _ 


^  +  2  J  c„  Ri  (u  -u  ) 

2  g  c  e' 


(i+0) 


Recognizing  that 


and 


u  -u  =  (l-e‘^/^®)  (u.-u  ) 
c  e  /  ^  j  e'' 


il  u 

s  ^ 


2  t'  dt' 


(or,  if  n  =  ;lf/(|r 

J 


u 

P  e 

2  ~  n  dn  )  , 
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ve  way  write  (4o)  as 


X  = 


2 


a[u^/(u.-u^)]  ds' 


Co(l-e 


ni  (s') 

-is'-  -1'^  ^ 


)  ( 

u 


o  e  '  j 

2  —  n  dn  ) 
u  ' 


(hi) 


WoWj  from  (35 ), 


thus 


U  -  +  (u^-u^)(l-e‘®®)  P*  (s,n); 


<as' 


ni (s' ) 

1q'  n-t 


-is ' V  * 


Cg(l-e  )  J  [  (l-e  )  p"  +  u  /(u.-u  )]"■''  2n'  dn' 

^  J  ^ 


m 


Since  ilf 

J 


i 


O^P  u  ''  ’ 
e  e 


u.  i  p  i 


X  =  X  + 


-  (5^)'  , 


o  2c  ^u.-u  '  V 


J  e  e 


^  u.-u 
J  e 


(^3) 


where  F 


ni  (s') 

'"2 


lo'  2 


2n '  dn ' 


■ts'^  % 


(1-e-^^  )  P 


u .  -u 
J  e 


dS 


The  evaluation  of  the  integral  becomes  somewhat  easier  far  from  the  origin^,  where 


Ik 


and 


Then 


P*(S,n)  cs  6"“^  . 


S  in2  ' 

or  p 


2 

2S 

2n'dn' 

Q) 

[U 

4 

CO 

1 

0) 

1 — 1  1 

2S' 

(44) 


u 


If  we  let  k  =  -  ,  the  inner  Integral  may  be  evaluated  analytically  to  give 


u  .-u 
J  e 


F  =  F(S^)  + 


r^_  i 

^2ks'  ^ 


2k  +  ^  4 


In 


2S 


k  + 


■}] 


dS 


(45) 


2S 


The  radial  transform  to  the  physical  plane  is  given  as 

f  p  u 

r  =  [2  J  (46) 

o 

Because  In  an  afterburning  rocket  exhaust  the  density  is  not  a  simple  function  of 
the  velocity,  this  transform  requires  a  greater  knowledge  of  the  details  of  the  flow 
field.  We  shall  discuss  the  requirements  in  a  later  section. 


C^se  A  (u  =  O): 

e  ' 

In  this  case  we  express  the  von  Mises  transform  as 

pur  =  p  .Uj  ♦  pvr  ,  -  Pju.  ♦  g  , 


(ij-7) 
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and  the  momentum  equation  becomes 


du  1  5|-/^2  2/  2  ,2»,  Su-i 


m 


With  the  same  assumptions  about  the  constancy  of  s_^  and  the  reference  density 


and  the  same  sort  of  linearization^  (i+8)  becomes 


5u  5  r .  Su-, 

3x  "  p  .  u .  i|;  3;j;  ^  SijrJ 
J  J 


m 


By  introducing  a  new  axial  coordinate 

?  - 


X  -1  * 

P  A  e 

u  .r . 

J  <] 


(50) 


the  heat  conduction  equation  is  obtained  as  before. 


^  _  J.  A 


(51) 


In  this  case  the  solution  is 


u 

u . 


A  ,  1 


=  (1-e  'J  )  P  ( 


M/r. 

J  J 


(52) 


To  transform  the  solution  to  the  physical  plane  we  again  assume  that  in  the 
upstream  region  of  the  jet 


s  =  c  ^  X  u  .  /ct 

V  1  j' 


and  downstream 


e  =  c^  Ri  u  /a 
V  2  t  c/ 


Then  in  the  upstream  region  we  obtain  the  solution 


l6 


(53) 


J  ""j  ""l 

which  my  he  recognized  as  the  limiting  case  of  Equation.  (38)  as  -•  0, 
In  the  downstream  region  we  have 


X  =  X  + 

o 


p  .  CTU.r  .d| 
J  J  J 


(5k) 


Proceeding  as  before,  we  find  that 


or  p  ^ 
o  2c2 


dS 


r 

nni(S)  -1 

u 

1q  ' 

2  /  / 

2n  dn 

P  (S,n) 

(55) 


It  will  be  noted  that  (55)  is  the  limiting  case  of  (h3)  when  u  ->  0.  Thus 

e  ■* 

or  p  ^ 

^  ^  ^p  )  F(S^,S,0)  .  (56) 

Far  from  the  origin,  (56)  may  be  shown  to  become 

or  p  I 

^  =  ^o  +  ^  [F(S^,S^,0)  +  S-S^]  ,  (57) 

which  may  be  derived  from  (56)  by  introducing  the  limiting  form  of  P*  or  from  (l<-5) 
by  taking  the  limit  as  k  -*  0.  The  radial  transform  is  performed  as  before. 

3.  DESCRIPTION  OP  FLOW  CHEMISTRY 

In  a  reacting  turbulent  flow,  a  complete  description  of  the  progress  of  chemical 
reactions  would  be  very  complex.  The  rate  expressions  describing  mean  reaction 
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velocities  'become  (like  other  expressions)  averages  over  fluctuating  quantities^  and 
the  proper  method  of  taking  this  average  so  as  to  produce  correct  results  has  not 
been  discovered.  We  therefore  have  essentially  three  alternatives :  the  reaction 
effects  may  be  neglected;  the  reactions  may  be  accounted  for  by  using  finite-rate 
expressions  ■which  are  adjusted  to  give  good  agreement  "With  experimental  results;  or 
the  reactions  may  be  assumed  to  proceed  very  rapidly  so  that  a  condition  of  Ic'al 
chemical  equilibrium  exists .  The  first  alternative  cannot  account  for  the  large 
energy  release  which  occurs  during  afterburning  of  the  exhaust  gases.  The  second 
requires  extensive  experimental  data  covering  many  different  conditions  in  order 
that  the  "phenomenological"  expressions  may  be  checked.  The  third  has  drawbacks  of 
its  own^  which  will  be  discussed  below ^  but  still  seems  the  most  reasonable  procedure 
for  a  rocket  exhaust  at  low  altitude.  It  is  this  approach  that  is  used  here. 

Even  if  reactions  are  very  fast,  it  is  strictly  improper  to  employ  mean-flow 
variables  in  the  equilibrium  expressions  unless  all  the  fluctuations  about  the  mean 
values  are  small.  However,  one  may  hypothesize  that  some  distance  downstream  from 
the  primary  reaction  zone  the  heat  released  by  combustion  of  a  given  element  of  fluid 
will  have  been  spread  by  diffusive  effects  to  neighboring  elements,  so  that  the  mean 
density  in  the  region  considered  is  approximately  that  given  by  assuming  that  the 
fluid  is  locally  "well-stirred".  We  therefore  should  expect  that  the  local  equil¬ 
ibrium  solution  described  by  conventional  thermodynamic  expressions  in  terms  of  mean 
flow  variables  represents  an  approximation  which  is  increasingly  good  as  one  moves 
downstream,  away  from  the  region  where  exhaust  gases  and  the  atmosphere  first  meet. 

We  have  treated  the  chemistry  of  the  reacting  plume  after  the  "adiabatic  flame 
temperature"  approach  proposed  by  Boynton  and  Heu.^^  Here  the  density,  composition, 
velocity,  and  temperature  may  be  expressed  as  functions  of  one  variable,  the  mass 
fraction  of  exhaust  gases  in  the  local  flow.  The  computation  assumes  local 
conservation  of  momentum,  mass,  and  total  enthalpy  (including  heats  of  formation)  on 
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mixing.  This  represents  a  solution  of  the  set  of  equations  (2)  if  the  fluxes  of 
mass  and  heat  and  the  shear  stresses  are  represented  by  equations  of  the  form 

T 


r,i  Sr 


T 


ST 


(78) 


T  au 

T  =  Pe  5-  i 

x,r  V  or  J 

and  if  in  addition  everywhere  in  the  flow  field.  Under  these  conditions. 


if  the  local  mss  fraction  of  fluid  which  originally  came  from  the  exhaust  jet  is  c., 

J 


then 


u  =  c  .u .  +  (l-c  . )  u  , 
JO  J  e 


h  +  i  u^  =  c  (h  +  (l-c  )(h  . 

J  J  J  J  “  ^ 


(79) 


The  method  of  Boynton  and  Neu  is  to  assume  that  u.  and  h.  are  determined  by  an 

0  J 

Isentropic  expansion  of  the  exhaust  gases  to  ambient  pressure.  (This  is  not  an 

essential  feature  of  the  computation,  and  a  non-isentropic  expansion  -  for  example, 

by  Thomson's  method  -  may  be  employed  if  the  nozzle  exit  and  free  stream  pressures 

are  greatly  different.)  For  a  given  value  of  c.,  the  values  of  u  and  h  are  computed 

J 

from  (79)*  From  the  thermodynamic  properties  of  the  various  chemical  constituents, 
the  temperature  of  the  mixture  is  found  by  a  straightforward  iterative  procedure .  At 
this  point,  no  chemical  reaction  is  assumed. 

The  effects  of  chemical  reactions  are  considered  in  the  following  manner,  A 
trial  temperature  is  chosen,  and  from  this  temperature  and  the  known  pressure  and 

atomic  composition  of  the  mixture  the  mole  numbers  (per  mole  of  exhaust  fluid)  are 

lb 

calculated  by  the  technique  of  free  energy  minimization.  This  process  is  then 
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Iterated  until  the  static  enthalpy  of  the  "burned  gas  at  the  trial  temperature  is 
equal  to  the  value  of  h  calculated  from  (79)*  It  is  assumed  that  the  energy  release 
does  not  accelerate  the  flow,  so  that  the  velocity  is  the  same  regardless  of  whether 
the  fluid  has  burned.  This  assumption  is  compatible  with  the  coordinate  transfor¬ 
mations  (l)  and  (27)  in  the  above  analysis  of  the  momentum  equation,  where  it  is 
assumed  that  the  only  effect  of  a  density  change  is  a  squeezing  or  stretching  of  the 
radial  coordinate .  Since 


u-u 

e 


(80) 


this  analysis  can  be  combined  directly  with  Libby’s  jet  flow  theory  to  provide  a 
map  of  velocity,  temperature,  and  composition  throughout  the  flow  field. 


4.  THE  COMPUTER  PROGRAMS  "SHARP"  AliH)  "SHEAR" 

The  analysis  presented  in  the  preceding  sections  has  been  embodied  in  two 
computer  programs.  The  first  program,  "SHARP",  computes  the  detailed  composition 
and  temperature  at  a  given  mixing  ratio  between  two  streams  whose  properties  are 
known,  using  the  techniques  of  the  previous  section.  The  second  program,  "SHEAR", 
uses  the  output  data  of  "SHARP"  to  compute  the  flow  field  by  an  integral  transformation 
of  Libby's  solution  to  the  physical  plane. 

The  program  "SHARP"  is  built  up  of  a  number  of  subroutines,  some  of  which  were 
developed  for  earlier  programs  ("Polyphase  Chemical  Systems,"  which  is  a  thermo¬ 
chemical  calculation  of  equilibrium  temperature  and  properties,  and  "PEEP,"  a 
propellant  performance  program)  and  some  of  which  are  specific  to  the  afterburning 
problem.  As  a  result,  much  of  the  input  information  required  by  SHARP  is  that  used 
by  PEEP.  The  following  is  a  listing  of  the  input  needed  for  SHARP  and  prescriptions 
for  obtaining  it . 

Card  No.  1  (Field  width  3):  NLS  -  number  of  chemical  species  in  the  mixed  flow; 

* 

A  condensed  phase  is  a  different’ species  for  this  purpose. 
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NT  -  number  of  temperatures  in  the  thermodynamic  property  tables j  HALT  -  the  number 
of  different  altitudes  to  be  considered;  NCAZR  -  the  number  of  atmospheric  constituents 
(usually  two,  and  Og).  One  card. 

ferd  Group  No.  2  (Field  width  l4.);  The  thermodynamic  property  tables,  one  set 
of  NT  cards  per  species.  These  contain:  T  -  temperature,  °K;  S  -  entropy  (third 
law  value)  in  cal/gm.  mole-°K;  H  -  enthalpy,  defined  as  the  heat  of  formation  at 
zero  K  (heats  of  foimation  of  the  elements  in  their  normally  occurring  states  - 
e.g.,  H2(g),  C(s)  -  are  taken  as  zero)  plus  the  heat  required  to  raise  the  species 
to  temperature  T,  in  cal/gm.  mole;  CP  -  heat  capacity  in  cal/gm.  mole-°K.  Additional 
Information,  such  as  species  identity,  may  be  included  at  the  right-hand  side  of  the 
card.  Cards  for  most  species  of  general  interest  are  available  in  the  Space  Science 
Laboratory.  NT  x  NDS  cards. 

&.rd  Group  No.  3  (Field  width  l^J-):  Atmospheric  and  trajectory  properties,  two 
cards  for  each  altitude.  First  card  contains:  ALT  -  altitude;  PALT  -  atmospheric 
pressure  in  atm.;  VALT  -  missile  velocity  in  ft/sec.;  TALT  -  atmospheric  temperature 
in  °K.  Second  card  contains  mole  fractions  of  atmospheric  constituents,  NCAIR 
entries.  Atmospheric  data  is  generally  obtained  from  standard  tables,  velocity 
data  from  a  trajectory.  2  x  NALT  cards. 

Card  Group  No.  4  (Field  width  3):  IMR  -  number  of  mass  ratios  of  air/jet  gases 
to  be  considered;  NDA  -  number  of  different  atoms,  or  of  some  other  ultimate 
components  of  the  system;  NLP  -  number  of  different  phases  in  the  system.  One  card. 

Ch,rd  Group  No.  5  (Field  width  l4):  TING  -  lowest  temperature  in  thermodynamic 
nroperty  tables;  TOE  -  initial  guess  to  temperature  of  expanded  gases;  TOB  -  initial 
guess  to  temperature  of  burned  gases;  HO  -  enthaloy  of  the  system  at  some  reference 
state  (generally  chamber  conditions)  in  cal/gm.  (This  may  be  obtained  by  taking  the 
enthalpy  of  one  mole  of  the  rocket's  fuel  in  cal/gm.  mole,  adding  to  it  the  enthalpy 
of  as  many  moles  of  oxidizer  as  are  required  by  the  mixture  ratio  at  which  the  rocket 
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IS  operated,  and  dividing  by  the  mass  of  the  fuel/oxidizer  system  in  atomic  weight 
units.).  VO  -  velocity  of  the  system  at  some  reference  state  in  ft/sec.  (zero  if 
chamber  or  reservoir  conditions  are  chosen).  One  card. 

Card  Group  No.  6  (Field  width  Ik):  FMR  -  the  mass  ratios  of  air  to  jet  gas  to 
be  considered.  IMR  numbers,  placed  5/card. 

Card  Group  No.  7  (Field  width  lU):  FMJ  -  the  atomic  weights  of  the  different 
atoms  in  the  system.  NBA  numbers  placed  5/card, 

&rd  Group  No,  8  (Field  width  3):  lEHASE  -  the  Identification  of  the  phase  in 
which  a  particular  species  helongs  (the  gas  phase  Is  always  phase  no.  1).  KDS 
number,  20/card. 

Card  Group  No.  9  (Field  width  6):  V  -  origlK,lly,  piolar  volumes  of  condensed 
species;  now  replaced  by  zeros,  NDS  number,  lo/card. 

Card  Group  No.  10  (Field  width  6):  A  -  number  of  atoms  of  the  jth  kind  in  the 


ith  molecule.  (For  example,  a 


'H,H20 


=  2;  a. 


=  1,  etc.)  NDS  entries,  lO/card, 


for  each  atom;  start  a  new  card  when  starting  a  new  atom.  MBA  groups  of  NDS  entries 
total. 

Card  Group  No.  11  (Field  width  ih) :  ZIN  -  composition  of  jet  gas  in  moles/ 
mole  of  "fuel”.  Mole  numbers  of  atmospheric  constituents  not  present  in  the  initial 
Jet,  or  of  species  produced  by  the  combustion  process,  must  not  be  zero,  but  are  set 
equal  to  a  very  small  number.  NDS  numbers,  5/card. 

Ckrd  Group  No,  12  (Field  width  l4) :  FM  -  molecular  weights  of  species.  NDS 
numbers,  5/card. 

Curd  Group  No.  13  (Field  width  6):  VAR  -  identlficntlon  of  the  epeclee  (Ng,  E^, 
etc.),  NDS  entries,  12/card. 

Card  Group  No.  l4  (ffey  use  whole  card):  NAME  -  identification  of  the  run.  One 

card. 

Ckrd  Group  No,  15  (Field  width  l4):  TO  -  temperature  (in  °K)  from  which 
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isentropic  exraansion  to  ambient  pressure  begins;  PO  -  pressure  (in  atm)  from  which 
expansion  begins.  These  may  be  chosen  more  or  less  arbitrarily  (i.e.,  rocket 

chamber,  throat,  or  exit  conditions),  but  must  be  consistent  with  each  other.  One 
card. 

Wherever  entries  occur  which  are  associated  with  specific  species,  or  specific 
atoms,  the  entries  must  always  be  made  in  the  same  order. 

The  output  from  SHARP  consists  of  the  following  quantities:  Altitude, 
atmospheric  pressure  PA  and  missile  velocity,  VA;  temperature,  TE;  and  velocity,  VE 
of  the  jet  gases  expanded  to  ambient  pressure;  mixing  ratio,  R;  temperatures  of  the 
unreacted  (tj)  and  equilibrium  (TB)  jet/air  mixture  and  mixture  velocity  VJ;  mole 
numbers  of  the  species  on  mixing  and  after  reaction.  This  information  is  repeated 
for  each  mixing  ratio  and  each  altitude. 

The  program  "SHEAR”  is  the  second  version  of  a  jet  mixing  program.  An  earlier 
version,  developed  for  an  application  to  a  launch  phase  countermeasures  study, 
contained  a  logical  error  which  affected  the  axial  coordinate.  In  addition,  the 
program  was  quite  inefficient  because  a  constant  step  size  was  used.  In  the  present 
version  the  step  size  of  both  the  radial  and  axial  transformation  is  controlled  by 
means  of  a  check  against  the  properties  of  the  solution  in  the  (§,Y)  plane  and  an 
interval  specified  by  the  user.  This  results  in  a  variable  spacing  in  the  output. 

The  input  required  for  SHEAR  consists  of  the  following  Information: 

ferd  Group  Wo.  1  (Format  2112,  lfE12.1+,  1E8.2);  WDS  -  number  Of  species;  WC  - 
number  of  Jet/air  ratios  tabulated;  TEST  -  the  value  of  the  jet  mass  fraction  at 
which  the  transition  between  the  "upstream"  and  the  "downstream"  eddy  viscosity 
model  is  made;  TOL  -  the  fractional  amount  by  which  the  centerline  mass  fraction  is 
diminished  to  define  the  axial  step  size  A§;  PRINT  -  a  number  which  is  set  negative 
for  no  printout  of  input,  zero  for  printout  of  Input  except  for  P*  tables,  and 
positive  for  printout  of  all  input;  ENDR  -  the  air  mass  fraction  at  which  the  radial 
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transformation  is  terminated;  ENDX  -  the  centerline  jet  mass  fraction  at  which  the 
axial  transformation  is  terminated.  One  card. 

Card  Group  No.  2  (Field  width  12):  "  value  of  jet  mass  fraction;  T.  -  value 

of  temperature  for  the  particular  value  of  R;  Z. .  -  values  of  species  mole  numbers 
(in  any  units,  moles/mole  of  mixture,  moles/mole  of  fuel,  etc.)  for  the  particular 
value  of  R,  in  consistent  order.  6  numbers/card  until  all  NDS  species  are  accounted 
for,  therefore  N  cards/value  of  R.  Start  a  new  card  for  the  next  value  of  F.;  NC  x 
N  cards  total. 

Card  Group  No.  3  (Field  width  12):  M  -  molecular  weights  of  the  species. 

NFS  entries,  6/card. 

Card  Group  No.  4  (Field  width  12);  P  -  ambient  pressure;  TP  -  temperature  of 
jet;  TE  -  temperature  of  air  or  free  stream;  A  -  radius  of  jet;  UP  -  jet  velocity; 

UE  -  air  velocity;  SIGMA  -  Jfech  number  correction  to  eddy  viscosity;  Cl  -  upstream 
region  eddy  viscosity  constant;  C2  -  downstream  region  eddy  viscosity  constant.  Two 
cards.  Any  absolute  units  may  be  used  for  pressure  and  temperature,  and  any  units 
for  velocity  and  jet  radius. 

The  SHEAR  program  also  requires  as  input  the  values  of  the  P^*-  function.  A 
tabulation  versus  a  =  ,  and  P  =  (l-t/tj  )/V'2§/tj'  has  been  prepared  from  the  data 

given  by  >fe.sters  and  stored  on  tape  5274.  This  tape  is  available  in  the  Digital 
Computer  laboratory. 

The  output  of  SHEAR  consists  of  the  following  items;  N,  the  index  of  the  radial 

step  (N=1  is  always  the  centerline);  X,  the  axial  coordinate  in  the  same  units  used 

for  A;  R,  the  radial  coordinate  in  the  same  units  used  for  A;  T,  the  local  temperature; 

U,  the  local  velocity;  and  the  local  composition  in  mole  fractions. 

SHEAR  also  produces  punched  card  output  consisting  of  X,  R,  T,  U,  and  the  mole 

fractions  of  the  first  four  species.  This  information  is  presented  in  a  8E10.4 

format.  These  cards  may  be  used  as  input  to  any  calculations  to  be  made  using  the 
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flow  field  data.  One  computer  run  produces  about  JOO  cards,  depending  on  the 
tolerance  used. 

The  operating  time  of  these  programs  is  apparently  quite  fast.  SHARP  load 
times  average  5  millihours  on  the  7094;  execution  time  for  12  mixing  ratios  of  an 
RP-1  jet  with  the  atmosphere  at  one  altitude  is  J  millihours.  The  corresponding 
SHEAR  run  had  an  execution  time  of  32  millihours .  It  is  therefore  possible  to  solve 
a  complete  far-field  exhaust  flow  problem  in  less  than  50  millihours  of  machine  time 
once  the  input  data  are  known. 

The  analysis  of  the  turbulent  mixing  problem  is  equally  applicable  to  the 

turbulent  axially  symmetric  wake.  In  the  programming  of  SHEAR,  the  eddy  viscosity 

was  set  proportional  to  the  absolute  magnitude  of  the  velocity  difference  u-u  so 

e  ^ 

that  SHEAR  may  be  used  for  a  wake  problem  where  u^  may  represent,  for  example, 
conditions  at  the  "neck". 

5.  SUMMARY 

At  this  point  it  seems  desirable  to  summari2e  the  foregoing  analysis  and 
discussion.  We  are  Interested  in  describing  the  flow  field  of  a  rocket  exhaust  jet 
at  low  and  moderate  altitudes .  The  most  important  processes  governing  the  flow  are 
the  turbulent  mixing  of  the  exhaust  gases  with  the  atmosphere  and  the  combustion  of 
the  fuel-rich  exhaust.  The  method  of  analysis  follows  closely  the  analysis  of  Libby 
developed  for  the  supersonic  combustion  of  hydrogen.  The  value  of  the  eddy  viscosity 
appropriate  to  a  jet  flow  is  derived  from  an  analysis  of  the  conditions  under  which 
compressible  jets  exhibit  self -preserving  behavior.  It  is  found  that  the  spreading 
coefficient  is  independent  of  density  when  expressed  in  the  proper  coordinate  system, 
but  that  an  effect  of  Ifech  or  Crocco  number  appears  to  be  present,  (it  is  shown  that 
this  effect  can  be  related  to  the  loss  of  turbulence  energy  as  radiated  sound.)  A 
consistent  development  of  Libby's  theory  is  given,  allowing  smooth  transitions  between 


25 


the  ease  of  zero  and  finite  free-stream  velocity.  It  is  shown  that  the  development 
of  the  turbulent  jet  depends  on  the  density  ratio  p./p  ^  the  velocity  defect  ratio 

viscosity  constants  c^^  and  o.  The  effects  of  reaction 
between  the  exhaust  and  the  atmosphere  are  accounted  for  by  assuming  that  the  flow  is 
locally  well-stirred  and  in  chemical  equilibrium.  Effective  turbulent  Prandtl  and 
Schmidt  numbers  are  assumed  to  be  identically  one .  Under  these  assumptions  the 
chemical  problem  may  be  decoupled  from  the  hydrodynamical  problem. 

Two  computer  programs  have  been  prepared  which,  when  used  together,  solve  the 
afterburning  rocket  exhaust  flow  problem  for  given  initial  conditions.  Instructions 
are  given  for  operating  these  programs.  In  the  first  appendix,  the  results  of  Libby's 
linearized  theory  are  compared  in  detail  with  well-known  self -preserving  solutions 
and  shown  to  be  nearly  equivalent.  Values  of  the  eddy  viscosity  constants  c^^  and  Cg 
are  derived  from  carefully  obtained  experimental  data.  In  the  second  appendix,  a 
comparison  is  made  with  the  results  of  constant -density  jet  flow  experiments  and  with 
rocket  exhaust  experiments  conducted  in  a  wind  tunnel,  and  the  agreement  is  shown  to 
be  satisfactory. 
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APBENDIX  1 


COMPARISON  OF  SELF- PRESERVING  AM)  LINEARIZED  THEORY 


AM)  DETERMINATION  OP  AND 


The  effects  of  the  linearization  on  Libhy's  solution  can  be  examined  by 
comparing  the  predicted  axial  and  radial  velocity  profiles  with  those  developed  in 
self -preserving  (asymptotic)  Jet  flows.  Let  us  first  consider  Case  A,  u^  =  0. 


Along  the  axis, 


-r.H 

u  =  u.  (l-e  ) 


c  J 


(A-1) 


or,  as  I  becomes  very  large. 


(A-2) 


1 

CTr.  P  •  ^  2£ 

Now  from  (57)  >  have  x  =  x^^  +  (■^)  —  ^ 

2  ^e  ^J 


Cp  p  ® 

2  /  ev 


I  =  V  (X-X,)  , 

J 


(A-3) 


whence 


(A-4) 


This  is  very  nearly  (l2)  for  f(Tl)  ==  1;  it  is  (l2)  exactly  if  I?  =  Acp .  From  (20 ), 


we  have 


^1  =  J  ^ 


A-1 


Then  for  exact  correspondence  between  (6l)  and  (l2),  we  must  have  c  = 

2  6 

The  radial  solution  far  downstream  may  be  expressed  (in  the  asymptotic  limit 
of  the  P  f inaction)  as 


=  e 


(A-5) 


or,  with  the  use  of  (6l),  as 


u 


=  e 


~'l'  A  — (^)®(x-x,  )r 

X  j 


(A-6) 


The  transformed  (or  equivalent  Incompressible)  radial  coordinate  is  expressed 


as 


P,u, 


2  A-i  A  dir' 

J  p  u  ^ 


P  -u. 
J  J 

P  U 
e  c 


^  (t^)'Aj  (^)^x-x^)r. 


(A- 7) 


With  a  little  manipulation,  (A-7)  becomes 


2  ^ 


which  may  be  rearranged  to  give 


l6c^ 


(A-8) 


where  T|  _  ctR/ (x  x.).  In  terms  of  a,  (65)  may  be  written 


-=  [1  +  ^]-^ 
U.  OCV-^ 


(A-9) 


A-2 


The  radial  profiles  obtained  from  (20 )  and  from  (A-9)  are  shown  as  functions 
of  Tl/2vff  in  Figure  1.  It  is  seen  that  the  theory  of  Libby  predicts  profiles  which 
are  narrower  than  the  self -preserving  profiles  near  the  center  of  the  Jet,  and  wider 
near  the  edge.  Therefore,  we  should  expect  that  properties  of  the  flow  which  are 
very  sensitive  to  the  radial  profiles,  and  particularly  to  the  profiles  at  the  edge 
of  the  Jet,  are  somewhat  inaccurately  described  by  the  linearized  theory  when 
Ug  =  0.  When  is  finite,  and  particularly  when  u  »  u.-u  ,  Libby's  theory  should 
be  as  accurate  as  the  similar  solution,  since  here  the  similar  solution  is  likewise 
linearized. 

Here  we  have 


err .  u .  p  P.if 

X  =  X  +  ^  (-^) 

1  2c^  'u.-u  '  'p  ' 

2  j  e  ^e 


2kG'  ds" 

k  + 


(A-10) 


If  A  »  1  while  S  «  1,  then  (A-10)  becomes 


o_  (  J  ^  /  J 


X  =  x„  + 

2  u.-u 


,  p  V  3/2 
j  e  e 


(A-11) 


where  the  lower  limit  to  the  integral  has  been  Incorporated  in  x^ .  Therefore 


3c,  2/3  (u,-uj2/3  /3 


lif.  "  ^2/jl,n2  ‘  /  V 


x-X2^2/3 


173  (5?f) 
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(A- 12) 


and 


u  -u 
c  e 

u.-u 
J  e 


2/3 


P/,  2  (u.-u  p 

'  3c_  3  '■n  ^  ) 


V3  „  1/3  2/3 


2  (u.-u  ) 

3  e' 
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X-X^ 


(A-13) 


In  this  case,  (A-13)  reduces  to  the  self-preseinring  expression  (23)  if 

A-3 


(A-11+) 


(u,-u  )2/3  ,2/3  p  1/6  2/3 

1=2  J  e  J  /ex  .  -5^2  x 

2  „  V3  ^ 

a  J 


Now  consider  the  radial  profile.  In  this  case,  in  the  far  downstream  region 


u-u 


u  -u 
G  e 


=  e 


-rAPj 


(A-15) 


and 


= 


P  u 

^  e  e  ,  I  / 

2  — —  d\k 

p  u  ^ 
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(A-l6) 


Upon  substituting  (A-I5)  into  (A-16)  and  integrating,  we  obtain,  after  a  bit  of 
manipulation. 


R 


2 


Therefore 


u-u^  -r2/2|^ 

-  =  e  , 

u-u  > 

c  e 


and  substituting  for  ?  from  (A-12),  we  find 

2/3  ,  1/6 


2/3 


f(Tl)  =  exp[-A/2(u.-u  (!£)  '  1 

J  e'  ^^T/3  V/  J  ^ 


(A-17) 


(A-18) 


(A-19) 


where  we  have  employed  the  relations 


and 


11  ^ 


R. 

J 


r . 

J 


(pj/p,)*  . 


Then  we  have  from  (A-I9)  and  the  definition  of  I  , 


^  = 


f(ll)l]dTl  =  2  —-3 


(u,-u  )2/3„2/3  1/6  3  2/3 


u 


(A-20) 


which  is  in  agreement  with  the  requirement  (A-14-)  for  if  Libby's  solution  and  the 
similar  solution  are  to  give  the  same  axial  decay.  Rote  also  that  (A-I9)  is  of  the 
same  form  as  (24),  and  we  have  therefore  defined  the  constant  a  in  terms  of  known 
quantities . 

Since  we  know  the  relation  that  must  exist  between  c^  and  cy,  the  constant 
appearing  in  (2o),  in  order  that  Libby's  analysis  and  the  similarity  analysis  give 
the  same  axial  velocity  decay,  we  may  determine  Cg  from  the  many  experiments  which 
have  been  performed  on  fully  developed  incompressible  jets.  In  Hinze's  book,  the 
value  of  O'  was  given  as  O.OOI96,  from  which 

Cg  =  =  0.0177. 

This  value  of  Cg  assures  the  proper  axial  decay  far  downstream.  It  should  be  noted 
that  this  number  is  quite  close  to  the  factor  relating  e,  u,  and  the  displacement 
thickness  in  the  outer  regions  of  a  turbulent  boundary  layer. 

In  the  present  case  we  have  two  slightly  different  values  of  R^  in  the  self¬ 
preserving  solution  and  in  the  linearized  solution.  Letting  y  =  Tl/vW,  we  may 
solve  (20)  and  (A-9)  for  y^.  The  results  are 

y|  =  0.910 
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for  the  self -preserving  expression  (20),  and 


=  0,816 

for  the  linearized  expression  (A-9),  Since 


81  =  y|  (x-x^) 

=  2y|  (x-x^), 

where  x^  is  the  virtual  origin  of  the  jet,  and  since  x  is  related  to  S  =  2|/ili.  by 
(57),  as 

R. 

we  have 


Ri  =  yi  s  7^ . 

2  2 

•X" 

In  order  to  give  the  correct  value  of  at  a  given  axial  station,  the  value  of  Ri 
derived  from  Libby's  solution  should  be  multiplied  by  O.9IO/0.816,  or  l.ll4. 
Therefore  the  effective  value  of  is 

c^  =  0.0197 

for  the  case  of  a  stagnant  ambient  medium. 

As  an  anproximate  interpolation  formula  between  the  two  cases  of  =  0  and 
(u.-u  )  «  u  ,  we  may  use 

c^  =  0.0177  +  0.0020  e"^  , 

where  k  =  u  / (u  -u  )  and  is  zero  when  u  =  0  and  infinite  when  u.-u  ^0.  Use  of  this 
e  j  e  e  j  0 

formula  will  give  the  proper  velocity  decay.  Since  the  apparent  turbulent  Schmidt 

number  found  in  incompressible  flow  is  approximately  0.7,  Cg  (and  also  c^^)  should  be 
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multiplied  by  this  value  if  composition  decay  is  of  interest.  No  really  valid 
approximation  exists  for  temperature  decay  in  a  system  where  viscous  dissipation 
effects  are  large,  but  it  might  be  expected  to  decay  more  like  composition  than 
velocity. 

The  behavior  in  the  upstream  portion  of  the  jet  is  determined  by  the  choice  of 
Cj^  and  the  value  of  §  at  which  the  transition  is  made  between  the  two  models  for  e*. 

V 

The  value  of  c^  can  be  determined  from  experimental  studies  of  mixing  layer  flows; 
here  the  solution  for  an  incompressible  flow  is  given  by  Schlichting  as^'^ 

u  =  i(u  +u^)[l-  erf  (P  |)]  ,  (A.21) 

j  e 

c  u.-u  c, 

where  p  =  i  (_  and  ^  is  defined  as  in  (36). 

J  e 

The  solution  obtained  by  Libby,  in  the  limit  that  ^  becomes  very  small,  is 
given  by 


u  =  u^  +  i  \l-eTf\ 


_ J 

^/§?7f7 

0 


(A-22) 


Evaluating  the  incompressible  radius  R,  we  find  that  when  the  linearization  is  valid 
R  «  ij;.  Defining  R  as  ijr  (u  /u  )  or  as  r  p  /p  we  may  rewrite  (A-22)  as 

dcl'^cl  ' 


U 


U 

R-R.  — 


u 

£ 

J  u, 


R.  — 


(A-23) 


Now  we  recognize  that  from  (38) 


/c,  u.-u 

J  V  O  U  R 
J  J 


Substituting  this  expression  in  (A-23),  and  assuming  u  /u.  —  1,  we  obtain 

®  <] 
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(A-24) 


u 


(A-25) 


u.-u 

j  e 

u. 


■]} 


and  by  Identifying  y  with  (R-Rj),  noting  that  u .  „  i  (Uj+n^),  and  defining  B  as  In 
Schlichting’s  solution,  we  find 


^  ^  ^  (A-26) 

This  IS  exactly  equivalent  to  (A-2i),  as  may  he  seen  by  recognizing  that 

Ug  =  i  [(Uj+u^)-(u^-u^)]  . 

Since  the  solutions  of  Libby  and  Schlichting  can  here  be  made  exactly  equivalent 
the  value  of  obtained  from  fitting  mixing-layer  experiments  to  the  Schlichting 
solution  can  be  used  in  transforming  Libby’s  solution  to  the  physical  plane. 
Schlichting  gives  a  value  of  p  =  13.5  as  the  best  fit  to  experimental  data;  Vasiliu 

selected  a  best  value  of  12.0.  These  values  are  for  u^  =  0  at  low  jet  Ifech  numbers, 
so  that 


C-,  =  ~  =  O.OOI7A  . 

The  value  of  can  be  determined  either  by  numerical  experiments,  or  by  choosing 
* 

It  so  that  as  given  by  the  upstream  expression  is  equal  to  e*  as  given  by  the 
downstream  expression. 


Schlichting’s  solution  is  also  linearized,  since  it  represents  the  first  two  terms 


of  a  power  series  in  (u.-u  )/(u  +u  ) 

J  ^  j 
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COMPARISON  OF  EXACT  AND  LINEARIZED  RADIAL 
VELOCITY  PROFILES  (  k  »  0  ) 
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RADIAL  PARAMETER,  17  /4  a 


APPENDIX  2 


COMPARISON  WITH  EXPERIMENT 


The  analysis  given  in  this  report  has  been  checked  against  experimental  data 

presented  by  Corrsin  and  Uberoi,^^  by  Forstall  and  Shapiro^^^  and  in  the  final 

20 

report  of  the  Ballistic  Missile  Plume  Wind-Tunnel  Research  Program.  The  comparison 
with  Corrsin  and  Uberoi  is  shown  in  Figure  2;  except  for  a  slight  axial  displacement ^ 
the  agreement  is  excellent.  (This  result  should  have  been  expected,  since  the  value 
of  Cg  taken  from  the  preceding  section  is  taken  from  this  experimental  work.)  The 
data  given  by  Forstall  and  Shapiro  show  considerable  experimental  scatter.  The 
results  of  a  comparison  with  data  from  two  different  velocity  ratios  are  shown 

in  Figures  3  and  k.  The  agreement  appears  to  be  fair^  considering  the  scatter  in  the 
data . 


Comparisons  with  data  from  actual  rocket  firings  are  difficult  to  obtain,  since 
only  a  few  data  exist.  The  only  source  of  data  at  altitude  seems  to  be  the  wind- 
tunnel  experiment  performed  by  Convair  in  the  AEDC  facility.  For  a  comparison  here, 
we  have  chosen  the  results  of  the  composition  samplings.  A  comparison  with  the  axial 
decay  data  at  sea  level  may  be  found  by  plotting  the  reciprocal  of  the  mass  fraction 
of  Jet  fluid  versus  axial  distance.  The  self -preserving  theory,  and  that  of  Libby 
in  its  asymptotic  form,  predict  a  linear  relationship  which  should  be  the  same  for 
all  Jets  when  plotted  versus  an  axial  coordinate  in  units  of  equivalent  Jet  radii, 
where 


r 

equiv 


The  results  of  such  a  plot  are  shown  in  Figure  5^  where  a  linear  relation  is  shown 

A-10 


to  exist.  We  also  show  the  reciprocal  velocity  decay  of  Corrsin  and  Uberoi,  and  the 
results  of  a  "SHEAE  II”  run  with  the  value  of  multiplied  by  O.65,  which  is  close 
to  the  value  of  the  turbulent  Schmidt  number  found  in  constant -density  Jets.^  The 
agreement  in  slope  between  the  computer  calculation  and  the  experiment  is  excellent. 

The  experimental  results  show  an  apparent  offset  of  the  virtual  origin  of  the  flow. 

(This  offset  also  appears  in  the  runs  at  altitude,  and  may  reflect  an  error  in  locating 
the  sampling  probe.) 

Radial  profiles  calculated  and  measured  are  shown  in  Figure  6,  where  we  have 
compared  radial  "slices"  with  nearly  the  same  axial  mass  fraction  for  the  experiment 
and  the  calculation.  The  agreement  is  again  good,  with  the  observed  profiles 
slightly  wider  than  the  calculated  ones  because  of  the  Schmidt  number  effect. 

A  comparison  at  altitude,  where  the  free  stream  velocity  is  appreciable,  has 
been  made  using  the  experimental  results  for  the  runs  at  50,000  ft.,  with  an  external 
Mach  number  of  I.75.  These  data  are  from  two  different  engines,  one  with  an  area 
ratio  of  25/1  and  the  other  with  an  area  of  8/1,  but  with  throat  area,  chamber 
pressure,  and  mixture  ratio  equal.  The  balanced  jet  assumption  should  apply  to  these 
data,  and  indeed  the  axial  decay  profiles  shown  in  Figure  7  follow  the  same  curve 
within  the  experimental  errors.  This  decay  curve  can  be  fitted  by  the  analysis, 
where  a  is  based  on  the  velocity  difference  Ifech  number  defined  by  the  inner , 'stream 
properties  and  an  effective  Schmidt  number  of  O.65  is  employed.  Since  the  apparent 
offset  from  the  origin  is  found  in  these  data  also,  we  show  only  the  region  of  the 
jet  where  the  data  were  obtained.  This  is  in  the  region  of  the  Jnrgest  volume  of 
hot  gas  in  an  afterburning  rocket  exhaust,  and  the  calculation  is  therefore  suitable 
for  use  in  a  radiance  prediction.  The  agreement  of  calculated  and  experimental 
radial  profiles  is  shown  in  Figure  8,  and  appears  to  be  as  good  as  could  be  expected 
considering  the  scatter  of  the  data. 
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EXPERIMENTAL  AND  CALCULATED  AXIAL  PROFILES 
CORRSIN  AND  UBEROI  ISOTHERMAL  DATA  k=0 


A-I2 


AXIAL  DISTANCE  .  x  /  r  : 


EXPERIMENTAL  AND  CALCULATED  AXIAL  PROFILES 
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WEIGHT  PER  CENT  OF  EXHAUST  GAS 
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EXPERIMENTAL  AND  CALCULATED 
RADIAL  PROFILES. 

WIND-TUNNEL  ROCKET  DATA.  SEA  LEVEL 


Figure 
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EXPERIMENTAL  AND  CALCULATED  AXIAL  PROFILES 
WIND-TUNNEL  ROCKET  DATA  ,  50,000  FT. 
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AXIAL  DISTANCE 


WEIGHT  PER  CENT  OF  EXHAUST  GAS 


Figure  8 
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